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ABSTRACT 

An important question in solar physics is whether solar microflares, the smallest currently observable flare events 
in X-rays, possess the same energetic properties as large flares. Recent surveys have suggested that microflares may 
be less efficient particle accelerators than large flares, and hence contribute less non-thermal energy, which may 
have implications for coronal heating mechanisms. We therefore explore the energetic properties of microflares 
by combining EUV and X-ray measurements. We present forward-fitting differential emission measure (DEM) 
analysis of 10 microflares. The fitting is constrained by combining, for the first time, high-temperature Reuven 
Ramaty High Energy Solar Spectroscopic Imager (RHESSI) observations and flux data from the Solar Dynamics 
Observatory (SDO) Atmospheric Imaging Assembly (AIA). Two fitting models are tested for the DEM; a Gaussian 
distribution and a uniform DEM profile. A Gaussian fit proved unable to explain the observations for any of 
the studied microflares. However, 8 of 10 events studied were reasonably fit by a uniform DEM profile. Hence 
microflare plasma can be considered to be significantly multi- thermal, and may not be significantly peaked or contain 
resolvable fine structure, within the uncertainties of the observational instruments. The thermal and non-thermal 
energy is estimated for each microflare, comparing the energy budget with an isothermal plasma assumption. From 
the multi-thermal fits the minimum non-thermal energy content was found to average approximately 30% of the 
estimated thermal energy. By comparison, under an isothermal model the non-thermal and thermal energy estimates 
were generally comparable. Hence, multi-thermal plasma is an important consideration for solar microflares that 
substantially alters their thermal and non-thermal energy content. 
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1. INTRODUCTION 

The size distribution of solar flares is well known to obey 
a power-law form (e.g., Drake 1971; Aschwanden & Freeland 
2012); thus, large powerful events occur much less frequently 
than small energy releases. However, despite this relationship, 
it is unclear whether small events exhibit energetic properties 
consistent with large flares. Solar microflares (see Hannah et al. 
2011 for a review), generally defined as events of GOES- class 
C, are the smallest currently observable events in the X-ray 
regime, and allow us to address this question. Recent surveys of 
microflares have suggested that, in many respects, microflares 
are similar to larger flares (Christe et al. 2008; Hannah et al. 
2008a, 2011). They occur in active regions and their frequency 
distribution is similar to large flares. Yet in one important 
aspect they differ significantly; their spectra above 10 keV 
(usually interpreted as non-thermal emission) are generally 
steep compared to large flares; interpreted as a power-law the 
spectral index of microflares is generally between —5 and —8 
(see Benz & Grigis 2002; Krucker et al. 2002; Christe et al. 2008; 
Hannah et al. 2008a), while for large flares it typically ranges 
from —2.5 to —4, (Saint-Hilaire et al. 2008), although there are 
some notable exceptions (Hannah et al. 2008b; O’Flannagain 
et al. 2013). A recent full review of microflare properties may 
be found in Hannah et al. (2011). One possible interpretation 
of this finding is that microflares are not as efficient particle 
accelerators as large flares. An important consequence of the 
steepness of microflare spectra is that the total energy in 
non-thermal electrons is strongly dependent on the lower energy 
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cutoff of the non-thermal spectra which is poorly constrained 
due to thermal emission at low energies. 

Historically, the thermal emission from microflares has been 
interpreted as isothermal despite the fact that it has been 
shown (e.g., McTiernan et al. 1999) that microflare plasma 
may be multi-thermal. Modeling a multi-thermal plasma can 
provide a more accurate estimate of the thermal energy content. 
Additionally, it also provides a more accurate measure of the 
low-energy cutoff and therefore of the non-thermal energy 
content. These two facts combined suggest that better insight 
into microflare energetics and therefore the energy input into 
the solar corona may be gained by considering multi-thermal 
models of the microflare plasma. 

In this work, we analyze the energy content of ten microflares 
by combining observations in X-rays and EUV from the Solar 
Dynamics Observatory (SDO) Atmospheric Imaging Assembly 
(AIA) and Reuven Ramaty High Energy Solar Spectroscopic 
Imager (RHESSI). Through these combined observations the 
thermal emission is modeled with a multi-thermal differential 
emission measure (DEM) instead of an isothermal model. We 
present an analysis of the energy content of these microflares 
and compare the estimated energy budgets with those obtained 
under isothermal assumptions. 

2. INSTRUMENTS AND DATA SELECTION 

Launched in 2010, the SDO (Pesnell et al. 2012) has provided 
a new perspective on solar activity. The AIA (Lemen et al. 2012) 
on board SDO observes the full Sun at seven EUV wavelengths; 
94 A, 131 A, 171 A, 193 A, 211 A, 304 A, and335A, coveringa 
temperature range from < 1 to > 10 MK. For each of these EUV 
wavelengths, AIA takes images with a pixel size of 0.6 arcsec at 
a cadence of 12 s, an unprecedented combination of spatial and 
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Figure 1 . Temperature response functions for the seven EUV channels of AIA 
(colored lines), and an equivalent temperature response function for RHESSI 
(black line). The RHESSI curve is scaled to one AIA pixel. The AIA temperature 
response functions include the empirical chi antifix correction, and are 
obtained using the SolarSoftWare (SSW) procedure aia_get_response . pro. 
(A color version of this figure is available in the online journal.) 

temporal resolution for full-disk solar imaging. Each of AIA’s 
observing wavelengths is associated with a featured, multi- 
thermal temperature response function. The 171 A channel for 
example, has peak sensitivity to plasma with a temperature 
of approximately T = 0.6 MK, whereas the 94 A channel is 
sensitive to plasma at both T — 1 MK and T = 6 MK, and the 
335 A channel is sensitive over a broad range from cool plasma 
as low as T = 100,000 K to hot plasma of T = 5 MK. The 
response functions are shown in Figure 1 . 

The most sensitive solar X-ray observations currently avail- 
able are provided by RHESSI (Lin et al. 2002). Launched in 
2002, RHESSI has provided high spatial and energy resolution 
observations of solar flares for over a decade, observing from 
3 keV-17 MeV. RHESSI has sub-second temporal resolution, 
and produces images of X-ray emission via rotating collima- 
tors, with a spatial resolution down to 2.3 arcsec. 

RHESSI provides excellent complementary data for temper- 
ature analysis of high-temperature solar plasma. As Figure 1 
shows, with the exception of 193 A the AIA channels are not 
very sensitive to plasma beyond T = 10-20 MK compared to 
lower temperatures, while RHESSF s lower sensitivity limit is 
T ~ 8 MK, and extends far higher than AIA (Figure 1, black 
line). Thus, RHESSI is ideal for constraining the properties of 
high-temperature plasma. In the energy range of interest for 
microflares (3-30 keV) and active regions, the energy resolu- 
tion of RHESSI spectral data is around 0.33 keV. 

In this paper, we focus on the smallest currently observ- 
able events with characterizable X-ray emission, so-called 
microflares. X-rays provide important information about the 
high-temperature plasma and provide direct observations of ac- 
celerated electrons present in microflares, allowing us to better 
constrain their energy content. For this work, we analyze a 
set of 10 microflares randomly chosen from the RHESSI flare 
catalogue. All of these events were observed by both SDO/ 
AIA and RHESSI between 2011 and 2012. To avoid saturation 
effects in AIA all of these events are of GOES- class B or lower. 
For reference throughout this paper, these events are numbered 
sequentially according to their date, with #1 representing the 
earliest event. Details of the event dates, times, and the chosen 
fit intervals are listed in Table 1 . 


3. ESTIMATING THE DIFFERENTIAL 
EMISSION MEASURE 

Several techniques are currently available to determine the 
DEM of hot solar plasma. The core challenge of determining 
a DEM is that it is an inverse problem which generally lacks a 
unique solution (see, for example, Ryan et al. 2014). Therefore, 
constraining the solution space using multiple instruments and 
physical understanding is important (see, e.g., Caspi et al. 2014). 

Most work with DEM distributions to date has been done 
for non-transient active region sources, using various methods, 
including regularization and maximum entropy, (Fludra & 
Sylwester 1986; Hannah & Kontar 2012), fitting with cubic 
spline functions (Brosius et al. 1996), Markov chain Monte 
Carlo (MCMC; Kashyap & Drake 1998), and forward fitting 
of single and multi- Gaussian functions (e.g., Aschwanden & 
Acton 2001; Aschwanden & Nightingale 2005; Aschwanden 
et al. 2013). DEM calculations for flares are less common, but 
can also be found, using forward-fitting and MEM (McTieman 
et al. 1999), regularization (Prato et al. 2006; Fletcher et al. 
2013), maximum likelihood (Kepa et al. 2008), and Monte Carlo 
methods (e.g., Reale et al. 2001) 

In this paper, we chose to apply a forward-fitting approach 
based upon that used by Aschwanden et al. (2013) on AIA 
images in the six optically thin wavelengths; 94 A, 13 1 A, 17 1 A, 
193 A, 211 A and 335 A. However, in this work we expand the 
fitting procedure to consider additional functions and also to 
simultaneously fit RHESSI spectra (see Figure 2). 

To estimate the DEM, fluxes for each wavelength are obtained 
from co-temporal AIA images. In order to study the same plasma 
with both AIA and RHESSI , first RHESSI images are generated 
at the event peak X-ray emission in the 6-12 keV channel 
using the CLEAN algorithm (Hurford et al. 2002). AIA data 
is averaged over the image integration time for the RHESSI 
image as well as over the area within the 50% contour of the 
RHESSI flare image. These intervals are listed in Table 1. Sun- 
integrated hard X-ray spectra for that same period were also 
generated, with appropriate background subtraction. 

As each AIA wavelength is sensitive to plasma at differ- 
ent temperatures (see Figure 1), each provides an additional 
data point for use in the DEM forward fit. For a given model 
DEM, the flux expected in each AIA wavelength may be cal- 
culated by multiplying this DEM distribution by the instrument 
response function for the given wavelength. Here we use version 
2 of the AIA temperature response function obtained from the 
SolarSoftWare (SSW) routine aia_get_response .pro, in- 
cluding the chiantif ix and evenorm correction factors. The 
expected flux is then compared with the real observed flux from 
AIA. The best fit is found by searching the parameter space (e.g., 
for a Gaussian DEM model peak temperature T p and Gaussian 
width cr) and finding the combination of parameters which pro- 
duces most accurately the flux values in the observed images. 
Simultaneously, the goodness of fit for each iteration of model 
parameters to the RHESSI count spectrum is determined as de- 
scribed in Section 3.1. 

In this paper, we explore two DEM models. First, we 
investigate the Gaussian model suggested by Aschwanden et al. 
(2013), which is characterized by a peak temperature T p and a 
Gaussian width factor cr, and can be written: 

/ (log T - log T„) 2 \ 

DEM = DEM 0 exp ( — V 6 6 p - ). (1) 

where DEMq is the peak emission measure of the distribution. 
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Figure 2. Schematic of the joint differential emission measure fitting procedure. AIA flux data from the six optically thin EUV channels is compared to the reconstructed 
flux from a model DEM function folded through the AIA temperature response function. Simultaneously, the RHESSI count spectrum is compared to a model count 
spectrum obtained by calculating the expected photon spectrum that results from the model DEM function, and utilizing the RHESSI response function to convert 
from a photon spectrum to a count spectrum. 

(A color version of this figure is available in the online journal.) 


Table 1 

Summary of Best-fit Results for Selected Microflares 


No. 

Event Date 

GOES 

Class 

GOES 

Peak Time 
(UT) 

RHESSI 

Fit Time 
(UT) 

Xc 

Gaussian 

[log T, a] 

Xc 

Uniform 
[log T, or] 

[7min, 7m ax ] 
(MK) 

1 

2011 Jun 5 

B3.5 

02:14 

02:13-02:14 

4.9 

6.60, 0.23 

3.2 

6.5, 0.67 

0.7, 14.8 

2 

2011 Jun 6 

B6.7 

13:20 

13:19-13:20 

7.2 

6.65, 0.21 

1.4 

6.6, 0.59 

1.0, 15.4 

3 

2011 Jun 21 

B2.8 

18:22 

18:22-18:23 

9.7 

6.55, 0.21 

1.1 

6.35, 0.73 

0.4, 9.5 

4 

2011 Jul 16 

B6.2 

17:04 

17:02-17:03 

5.9 

6.55, 0.27 

6.5 

6.45, 0.74 

0.5, 15.4 

5 

2011 Aug 26 

B4.2 

20:53 

20:53-20:54 

13.3 

6.60, 0.24 

1.9 

6.60, 0.59 

1.0, 15.4 

6 

2011 Oct 11 

B5.8 

00:35 

00:35-00:36 

19.4 

6.65,0.19 

2.0 

6.65, 0.46 

1.5, 12.9 

7 

2012 Jun 20 

B7.8 

15:53 

15:48-15:49 

12.8 

6.60, 0.24 

1.1 

6.6, 0.58 

1.0, 15.0 

8 

2012 Sep 10 

B8.7 

07:24 

07:22-07:23 

30.1 

6.60, 0.23 

7.2 

6.60, 0.56 

1.1, 14.5 

9 

2012 Sep 15 

B2.9 

22:44 

22:44-22:45 

24.5 

6.55, 0.23 

3.0 

6.50, 0.62 

0.8, 13.2 

10 

2012 Sep 27 

B7.5 

06:57 

06:56-06:57 

30.1 

6.60, 0.20 

3.7 

6.60, 0.50 

1.3, 12.6 


Second, to model a uniform emission measure we choose an 
Epstein function (see Nakariakov & Roberts 1995; Pascoe et al. 
2007; Inglis et al. 2009, for previous examples of use of this 
function), which may be written 


DEM = DEMosech 2 


( 


log T - log Tp 
a 



( 2 ) 


where a and DEM 0 are width and amplitude parameters as 
before, T p remains a temperature parameter denoting the center 
of the distribution, and n is a steepness parameter. For both 
models, the parameter search space consists of 40 T values 
equally distributed in log space between log T — 5.5 and log T = 
7.5, and 15a values equally spaced in the range 0.05-0.8. At 
each T and a, the amplitude parameter DEM 0 is estimated 
by minimizing the difference between the observed and the 
modeled AIA flux. 

This distribution is convenient as it allows a smooth variation 
between a classical Epstein profile, where n — 1, to a boxcar 
function where n oc (see Figure 3). Throughout this paper, 
we assume a boxcar approximation by setting n = 10. This 
choice tests a simple case for the functional form of the 
DEM, which is that an entirely uniform emission profile may 
adequately fit the combined AIA and RHESSI data, between 
low and high cut-off temperatures. The contention is not that 


such a simple DEM distribution represents the actual emission 
structure of these events, but that it serves as a useful first-order 
approximation. This hypothesis gives us insight into whether 
the observations require more complex forms of the DEM 
(e.g., Warren et al. 2012, 2013) and by extension additional 
parameters. Also, this function is able to fall off steeply at high 
T, an important consideration given previous studies of RHESSI 
spectra which have adequately fit high-temperature plasma with 
power-law or exponentially decaying emission (e.g., Hannah 
et al. 2008a). For convenience, we further define T max = T p + a 
and Tinin = T p — a , as these parameters provide a more intuitive 
description of a uniform emission profile than T p and a . 

Even for relatively small events such as microflares, saturation 
of AIA images can occur. Saturation is characterized by bleeding 
of flux from one pixel into neighboring pixels, the result being 
that the true flux information is spread over a substantial image 
portion. Significant image saturation can lead to inaccurate 
DEM fit results. To account for this, each microflare was 
examined during the chosen fitting interval to find the fraction of 
pixels under the 50% RHESSI contour in each AIA image where 
the data number exceeded 1.6 x 10 4 . This level corresponds 
approximately to the saturation level of AIA (Raftery et al. 
2011 ). 

Only two microflares showed any significant saturation, 
defined here as >10% of the relevant pixels. Event #4 (2011 
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Figure 3. Normalized Epstein function centered at x = 0, with steepness 
parameters of n = 1 (solid line), 3 (dashed line), and 10 (dot-dashed line). 
For comparison, a Gaussian centered at x = 0 is given by the red line. In all 
cases, width factor a — 1.5. 

(A color version of this figure is available in the online journal.) 


range is restricted to 4-8 keV since above 8 keV non-thermal ef- 
fects can be important. Additionally, for consistency, the DEM 
fits are performed on data from detector 1 of RHESSI in all 
cases. For additional simplicity, only the model spectrum de- 
rived directly from the DEM model is fit to the RHESSI count 
spectrum; other common corrections, such as albedo effects, and 
modification of the elemental abundances are neglected in order 
to reduce the number of free parameters, although it should 
be noted that including such small corrections would likely 
improve the RHESSI goodness-of-fit. To account for this, we 
allow for greater systematic uncertainty associated with the 
RHESSI data points, adopting a value of 4% throughout this 
work. This includes a standard 2% level uncertainty correspond- 
ing to uncertainty in the base RHESSI response function, com- 
bined with an additional 2% estimated uncertainty accounting 
for the small corrections listed above. 

Throughout this paper, we refer to the reduced x 2 values for 
AIA or RHESSI data alone defined by 


Xaia 


1 

Naia — n — 1 


X 

AIA 


{Jo - fm? 

a l IA 


(3) 


July 16) contained one saturated frame and event #10 (2012 
September 27) contained two saturated frames. These frames 
were removed from the analysis, meaning that the AIA flux 
data is averaged over fewer frames for these two events. 

Another issue is that the temperature response functions 
of SDO/AIA are associated with significant uncertainty that 
affects the best-fit DEM and subsequent energy estimates. 
These uncertainties remain the subject of active study (e.g., 
Aschwanden et al. 2013). In particular, the 94 A and 131 A 
channels are known to be affected at low temperatures by lines 
not currently accounted for in the chi ant i database. To account 
for this, we adopt a systematic uncertainty of 25% for the flux 
detected in each AIA channel as suggested by Boemer et al. 
(2012) and Guennou et al. (2012). 

3.1. Fitting RHESSI Spectra 


and 


Xhsi 


1 


A^hsi — n — 


tE 


(fo ~ fm ) 2 


HSI 


msi 


(4) 


where f Q is observed flux, f m is the expected model flux, a 
denotes the variance, n is the number of free parameters, and 
Na\a and Nu si represent the number of data points available 
for AIA and RHESSI , respectively. For both models used in 
this paper, the free parameters are log T p , a and EM {) (see 
Equations (1) and (2)), and hence the number of free parameters, 
n , is 3. 

The combined reduced x 2 is denoted by x 2 and is obtained 
as follows 


2 _ ~ f m f ( fo ~ fm) 2 

Xtotal — 2 Z Z2 + E2 

AIA ^AIA HSI ^HSI 


RHESSI spectra are critical in order to provide high- 
temperature constraints to the DEM model. Fits to the 
RHESSI count spectra were made using the OSPEX analy- 
sis package available through SSW. New fitting routines were 
developed for this investigation, multi_therm_gauss and 
multi_therm_epstein, and have been made freely available 
within OSPEX, where their functionality is described in detail. 2 
The fundamental feature of these routines is that, for a given 
set of DEM function parameters, the expected photon spectrum 
is generated by assuming bremsstrahlung continuum emission. 
Subsequently, the RHESSI detector response matrix (DRM) is 
used to predict the count flux that would be observed by RHESSI. 
This flux is compared with the observed flux in each energy bin 
within the fitting range, and ax 2 minimization is performed. 

Typically, for microflares, RHESSI detects significant emis- 
sion only in the range E ^ 15 keV. In this study, we fit 
background-subtracted RHESSI spectra, where the fitting time 
interval spans 60 s. This length of time integration is necessary 
for small events such as microflares in order to accumulate suf- 
ficient counts to generate an accurate spectrum. Table 1 lists the 
fit and background intervals of each event. The fitting energy 


2 OSPEX documentation is available online at 
http://hesperia.gsfc.nasa.gov/rhessi2/ 


and 


Ac N-n- 1’ 

where N = Naia + Nusi is the total number of data points. In 
this study, ZVhsi is always larger than Naia, which equals 6 in all 
cases, whereas over the 4-8 keV energy range ZVhsi =12. The 
consequence is that the RHESSI data contribute more strongly 
to x 2 . 

4. RESULTS 

We analyze the full set of events listed in Tables 1 and 2, in 
each case carrying out joint fitting of the AIA flux data and the 
RHESSI count spectrum as illustrated in Figure 2 and described 
in Section 3. Both the Gaussian and Epstein profiles are used, 
and the results of joint fitting using each model are displayed in 
Figures 4 and 5, respectively. 

Figure 4 shows the performance of the Gaussian model, where 
three panels are plotted for each microflare. The leftmost panel 
in each case shows two DEM models. The first is a reference 
best-fit where only the AIA flux data has been used (red line), 
and the second is a combined fit utilizing both AIA and RHESSI 
data (blue line). The middle panel shows the ratio between the 
combined best- fit modeled AIA flux and the actual observed flux 
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(g) 2012/06/20 


(h) 2012/09/10 




(i) 2012/09/15 o<ev) ,o 9 t ( k) (j) 2012/09/27 



Figure 4. Results of joint forward-fitting of a Gaussian DEM profile to AIA fluxes and RHESSI count spectra. The events are ordered sequentially according 
to the event date. For each sub-image, the leftmost panels show the combined best- fit DEM as a function of log T (blue line). For comparison, the DEM which best 
fits the AIA data alone (red line) is also shown. The middle panels show the ratio of the best-fit model AIA flux to the observed flux, while the rightmost panels show 
the comparison between the model count spectrum and the observed RHESSI spectrum. 

(A color version of this figure is available in the online journal.) 


Table 2 

Summary of Energy Estimates for Selected Microflares 


No. dL/dt tdm Trad Multi + Thick a Iso + Thick b 



(10 24 erg s” 1 ) 

(s) 

(10 27 erg) 

U t h 

(10 29 erg) 

Fnth 

(10 29 erg) 

E c 

(keV) 

Unth/Uth 

u t h 

(10 29 erg) 

Gith 

(10 29 erg) 

E c 

(keV) 

Fnth/ Gh 

1 

10.1 

180 

1.8 

0.49 

0.18 

14.3 

0.37 

0.45 

0.39 

11.8 

0.86 

2 

20.4 

420 

8.6 

1.5 

0.81 

12.3 

0.54 

1.5 

1.5 

11.1 

1.00 

3 

7.1 

300 

2.1 

0.34 

0.09 

11.2 

0.26 

0.36 

0.11 

11.4 

0.30 

4 








0.45 

0.19 

11.5 

0.42 

5 

5.4 

480 

2.6 

0.39 



<0.01 

0.42 

0.43 

11.8 

1.02 

6 

21.6 

300 

6.5 

0.72 

0.14 

15.2 

0.19 

0.71 

0.66 

9.4 

0.93 

7 

16.4 

900 

14.7 

1.3 

0.47 

14.3 

0.36 

1.34 

1.23 

11.8 

0.92 

8 








0.95 

0.35 

9.6 

0.37 

9 

8.2 

180 

1.5 

0.47 

0.14 

10.5 

0.30 

0.55 

0.07 

12.7 

0.13 

10 

17.0 

240 

4.0 

0.64 

0.07 

9.9 

0.11 

0.63 

0.32 

5.2 

0.51 


Notes. dL/dt and L ^mit are derived using Epstein DEM profile in all cases. 

a Thermal and non-thermal energy calculated by using the best-fit multi-thermal model and fitting the remaining high-energy emission with a thick target 
bremsstrahlung model. 

b Thermal and non-thermal energy calculated using a traditional isothermal-plus-thick target spectral model. 


in each channel. Finally, the right panel shows the comparison 
between the combined best-fit model count flux, and the actual 
count flux observed by RHESSI. 

The Xc values for each microflare (see Table 1) indicate that 
these events are not well-represented using a single Gaussian 


DEM model, with x«? ~ 5 for event #1 the best value obtained. 
Figure 4 reveals that the combined best-fit DEM is consistently 
shifted to favor a narrower distribution when compared to best- 
fits obtained using AIA data only. However, this combined fit 
tends to provide a poor reproduction of both the AIA fluxes 
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(a) 2011/06/05 




(c) 2011/06/21 


(d) 2011/07/16 




(g) 2012/06/20 


(h) 2012/09/10 





Figure 5. Results of joint forward-fitting of an Epstein DEM profile to AIA fluxes and RHESSI count spectra. The events are ordered sequentially according to the 
event date. For each sub-image, the leftmost panels show the combined best-fit DEM as a function of log T (blue). For comparison, the DEMs which best fit the 
AIA data alone (red line) is also shown. The middle panels show the ratio of the best-fit model AIA flux to the observed flux, while the rightmost panels show the 
comparison between the model count spectrum and the observed RHESSI spectrum. 

(A color version of this figure is available in the online journal.) 


and the RHESSI spectra. Hence, this model appears to be 
inappropriate for constraining both high- and low-temperature 
plasma simultaneously. This is discussed further in Section 4.1. 

The result of applying the Epstein model to the events is 
presented in Figure 5. Here we observe mixed results. For events 
#2, #3, #5, and #7, the model can be said to fully describe the 
data, with x c 2 < 2 . These values may be re-cast in terms of 
the probability value (or p- value), where values of p below a 
set threshold (e.g., 0.05 for a 95% significance level or 0.01 
for a 99% level) indicate that the model should be rejected as 
a full description of the data. The x<? values for these events 
correspond to probability values of p # 2 = 0.14, p # 3 = 0.35, 
p #5 =0.02 and p # 7 =0.35, respectively. 

Of the fully fit events, a visual inspection of the RHESSI 
spectra in Figure 5 reveals that microflares #3 and #5 are in fact 
fully fit by the multi-thermal model at energies above 8 keV. 
This implies that the non-thermal component of emission is 
either absent in these events or is sufficiently weak to be entirely 
masked by the thermal emission. For events #2 and #7, there is, 
by contrast, evidence of a break above 8 keV to a power-law 
distribution, implying the presence of significant non-thermal 
energy. 

For the remaining events, four have a Xc value between 2.0 
and 4.0 (events #1 , # 6 , #9, and #10). For these, the chosen model 
does not fully describe the data. However, Figure 5 shows that 


the form of the RHESSI spectrum is reasonably well reproduced 
for these events, as is the AIA flux at the 131 A, 193 A, 211 A, 
and 94 A wavelengths. The 171 A and 335 A emission for these 
events is the least well reconstructed. For example, in event 
#9, the predicted 335 A flux is ~ 1 .75 that of the observed flux, 
while in event #10 the 171 A flux is underestimated by a factor 
of 3. Since both of these channels have a substantial temperature 
response component below log T = 6.0, one explanation is that 
the model DEM function is not fully accounting for this lower 
temperature emission. 

The remaining two events, #4 and #8, measure x<? = 
6.5 and x«? = 7.2, respectively. Hence the simple uniform 
DEM distribution applied here is insufficient to describe these 
events. Since a single Gaussian DEM distribution also proved 
inadequate, the use of more complex DEM distributions in these 
cases is justified. 

4.1. A Case Study: 2011 June 21 

In order to explore these results in detail, an in-depth study of 
event #3 is presented here. This microflare originated from active 
region AR1 1236 and occurred at approximately 18:20 UT, as 
shown in Figure 6 . Figure 6 (b) contains both AIA lightcurves 
in each waveband and RHESSI lightcurves in a number of 
X-ray channels. The fitting interval used is denoted by the 
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(a) Microflare images at each AIA wavelength. A RHESSI image is 
overlayed as red contour on the 94 A panel. 




(b) AIA (top) and RHESSI (bottom) lightcurves of the 2011 June 21 
microflare in several channels 

Figure 6. Summary information for the 2011 June 21 microflare from 18:20 UT. (a) AIA images of the microflare at each EUV wavelength, with a reconstructed 
RHESSI image overlayed on the 94 A image as a red contour. The RHESSI image was constructed in the 6-12 keV range using the CLEAN algorithm and shows 
emission above the 50% level, (b) AIA lightcurves of the microflare region and immediate surrounds, and RHESSI lightcurve of the event displayed in a range of 
energy bands. 

(A color version of this figure is available in the online journal.) 
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vertical dashed lines. Figure 6(a) illustrates the spatial extent 
of the microflare in both AIA and RHESSI , with the RHESSI 
image shown as a contour on the 94 A panel from AIA. 

4.1.1. Fitting AIA Fluxes in Isolation 

Many previous studies of DEM have relied on a single 
instrument only (see, however, Caspi et al. 2014). To understand 
the contrast between this and joint analysis, we constrain the 
DEM for event #3 using data from AIA only. For the popular 
Gaussian DEM model, the results of this process are displayed 
in Figure 7. 

The top panel displays the best-fit DEM, while the middle 
panel compares the observed AIA flux data with the model 
flux in each channel derived from the best-fit Gaussian model. 
These panels show that the flux from five of the six channels 
is reproduced to within the error bars, while the remaining 
channel — 171 A — has a slightly larger discrepancy. In this 
case Xaia ~ 1-0, corresponding to a good fit. The Gaussian 
parameters are: logEMo = 21.58 cm -5 K -1 , o = 0.45, and 
log T p = 6.50, equivalent to ~3.2 MK. 

However, although the best-fit function is able to reproduce 
the observed fluxes in each AIA channel, in the bottom panel 
of Figure 7, the resulting RHESSI count spectrum is compared 
with the observed background- subtracted energy spectrum. A 
difference of approximately two orders of magnitude is evident. 

Hence, although an excellent fit to the AIA data, this DEM 
dramatically overestimates the emission at high temperatures. 
This is a limitation of the Gaussian function, which is unable 
to decrease sufficiently quickly at high T when a is large (e.g., 
a > 0.2), and is already too broad to correctly fit the low- 
temperature plasma. 

We similarly apply the Epstein profile to the AIA data for 
event #3, without utilizing RHESSI. The result of this fit is 
shown in Figure 8. 

Here the DEM function again provides an excellent fit to 
the AIA flux data, with Xaia ~ 0.7 and all six channel fluxes 
reproduced within the error bar. Moreover, the bottom panel 
in this case shows that the expected count flux in RHESSI (red 
line) is of the same order as the observed flux (blue, dashed line). 
However, the best-fit model remains unable to truly reproduce 
the RHESSI spectrum. 

The fact that both the Gaussian and Epstein models can be 
successfully fit to the AIA flux shows that AIA alone provides 
cannot distinguish between the models. Furthermore, the bottom 
panels of Figures 7 and 8 imply that, regardless of the chosen 
model, AIA data is insufficient to properly constrain emission 
at high temperatures. 

This is further confirmed by carrying out AIA-only fits on 
each of the sampled events. For each microflare, the RHESSI 
counts estimated from the Gaussian model exceed the actual 
observed emission by at least one order of magnitude, usually 
two, similar to that observed in Figure 7. The Epstein profile, by 
contrast, provides RHESSI count estimates of the same order 
as the observed spectrum 80% of the time, the exceptions 
being events #8 and #9. However, the Xhsi va l ues remain poor 
(generally in excess of 10), indicating that the AIA data is 
insufficient to determine the high-temperature cutoff correctly. 

4.1.2. Combining AIA and RHESSI in a Joint Fit 

For event #3 we examine the joint AIA and RHESSI fits 
for both the Gaussian and Epstein models. The results of the 
Gaussian model are shown in Figure 4(c), and the Epstein model 
in Figure 5(c). 
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Figure 7. Results of applying a best-fit Gaussian DEM function to the AIA flux 
data alone for the 2011 June 21 microflare. Top panel: the best fit model DEM. 
Middle panel: AIA flux comparison, where the actual observed flux is denoted 
by the diamonds and the model flux associated with the best-fit Gaussian DEM 
is denoted by the crosses. The ratio of the model flux to the observed flux is also 
shown for each AIA channel. Bottom panel: comparison of the model RHESSI 
count spectrum that results from the best-fit DEM (solid line) to the actual 
energy spectrum observed by RHESSI. The observed spectrum was obtained 
between 18:22 and 18:23 UT on 2011 June 21 (dashed line). 

(A color version of this figure is available in the online journal.) 


The result of the combined fitting of the Gaussian model is 
significantly different from the AIA-only fit. The model RHESSI 
spectrum is closer (albeit still a poor fit, with Xhsi = 12.6) to 
the observed count spectrum, as a result of the best-fit model 
being much narrower (cr = 0.21). This narrower Gaussian 
means that the contributions to the total emission reduce much 
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Figure 8. Results of applying a best-fit Epstein DEM function to the AIA flux 
data alone for the 2011 June 21 microflare. Top panel: the best-fit model DEM. 
Middle panel: AIA flux comparison, where the actual observed flux is denoted 
by the diamonds and the model flux associated with the best-fit DEM is denoted 
by the crosses. The ratio of the model flux to the observed flux is also shown for 
each AIA channel. Here Xaia = 0.7. Bottom panel: comparison of the model 
RHESSI count spectrum that results from the best-fit DEM (solid line) to the 
actual energy spectrum observed by RHESSI (dashed line). 

(A color version of this figure is available in the online journal.) 


more quickly in the high T range, thus avoiding the unphysical 
excess contributions discussed in Section 4.1.1. However, the 
consequence of this is that the narrow Gaussian is no longer able 
to accurately reproduce the observed AIA fluxes, as evidenced 
in the middle panel of Figure 4(c). Hence, the model fits neither 


the AIA or RHESSI data well, and the overall goodness of fit, 
Xc = 9.7, is unacceptable. 

The Epstein profile in contrast is able to fit the AIA and 
RHESSI data simultaneously, as Figure 5(c) shows. Five of 
the six AIA fluxes are reproduced within the error bars, and 
the RHESSI spectrum is well reproduced, with the result that 
Xc = 1-1. 

For event #3, a constant DEM distribution as a function of T , 
combined with a steep fall-off at a high temperature, provides 
a good fit to the joint data set, particularly considering that 
higher order corrections to the RHESSI spectrum were neglected 
(see Section 3.1). By comparison, joint fitting using a Gaussian 
model was unable to simultaneously fit both the AIA fluxes and 
the RHESSI spectrum with a single consistent function. Both 
Figure 4 and Table 1 show that this limitation of the Gaussian 
function was replicated for all of the microflares in the sample, 
with Xc > 4.0 in all cases, and generally much higher. Hence 
this model cannot be used to understand plasma emission over 
a wide temperature range. 

By contrast, the results achieved with the Epstein model 
consistently improve upon those obtained with the Gaussian 
model, with / 2 < 2.0 in four events. In four additional events, 
2.0 < Xc < 4.0. Hence, based on the RHESSI and AIA 
data, most of the studied microflares may be adequately fit by 
a uniform DEM profile with a steep high-temperature cutoff, 
implying that any finer structure in the DEM cannot be resolved 
by these instruments. 

5. MICROFLARE ENERGETICS 
5.1. Radiative Losses 

Given a well-constrained DEM for a plasma, it is possible 
to estimate the radiated energy loss rate as follows (e.g., 
Aschwanden 2005; Ryan et al. 2012). 

— = f DEM(r) x A (T) dT (7) 

dt J 

In this equation, DEM(T) is the amount of emission mul- 
tiplied by event area A and hence in units of cm -3 K -1 , and 
A (T) represents the radiative loss function, which has been in- 
vestigated by many authors (e.g., Rosner et al. 1978; Mewe 
& Gronenschild 1981; Cook et al. 1989; Martens et al. 2000; 
Klimchuk et al. 2008; Dudik et al. 201 1). 

An estimate of A (T) may be obtained from the chianti 
database after choosing appropriate coronal abundances. Al- 
though this function is only weakly dependent on density (e.g., 
White et al. 2005), the radiative losses vary by an order of mag- 
nitude as a function of temperature in the range 0.5 < T <5 MK. 

The total radiated energy can be estimated by the integral of 
the radiative loss rate dL/dt over the duration t dur of the event. 
However, for microflares it is difficult to resolve the behavior 
of the DEM distribution as a function of time, given the need to 
integrate the RHESSI counts for up to 60 s in order to obtain a 
viable spectrum. Instead, we estimate the total radiated energy 
as L rad ~ dL/dt x £ dur , where £ dur is the duration of the microflare 
estimated from the event lightcurves. 

Estimates of dL/dt and L rad are obtained for 8 of the 10 
events, with events #4 and #8 discarded due to poor fits. 
We find that the estimated radiative loss rates range between 
5 x 10 24 -2 x 10 25 erg s -1 , while the total radiative losses L rad 
range between 2 x 10 27 — 1 .5 x 10 28 erg, with a mean value of 
~5 x 10 27 erg. The estimates for each event may be found in 
Table 2. 
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This can be compared with estimates obtained from 
background-subtracted GOES data, where an isothermal plasma 
and coronal abundances are assumed (White et al. 2005). To es- 
timate the background for each event, a featureless time interval 
prior to the event was selected and taken to be the background 
level. To estimate L ra d , the estimated radiative losses dL/dt 
were integrated over the same duration as listed in Table 2. 
We find that the GOES energy estimates are of the same or- 
der of magnitude as those obtained from the fitted DEM dis- 
tribution. In general, the GOES-e stimated energies are higher 
than the estimates obtained from the multi-thermal DEMs, by a 
factor ~2. 


5.2. Thermal Energy 

For an isothermal plasma at temperature T, the thermal energy 
is given by (e.g., Hannah et al. 2008a; Emslie et al. 2012) 


U± = 3n e k B T fV, (8) 

where / is the filling factor, n e is the electron number density, 
V is the plasma volume, and k B is Boltzmann’s constant. Since 
n e = y/EM/fV, this may be rewritten as 


U th = 3£ b tVEM/V. (9) 


A recent statistical study by Hannah et al. (2008a) analyzed 
more than 25,000 microflares in this way, finding that their 
median thermal energy at the peak of the event was of the order 
10 28 erg, with a 2cr range of 10 26 -10 3 ° erg. In this study, the 
median volume V was estimated as 10 27 cm 3 with a 2 a range of 
5 x 10 25 -2 x 10 28 cm 3 , and the median emission measure was 
found to be 3 x 10 46 cm -3 with range 4 x 10 45 -2 x 10 47 cm -3 . 

The thermal energy may also be estimated for a multi-thermal 
plasma. In this case, it is necessary to account for all the energy 
at all T. From Equation (9) we may write 


f/th = 3 k B V 


1/2 f T DEM x TdT 
EM 1/2 


( 10 ) 


where DEM = dEM/dT , the DEM (in cm -3 K -1 ), EM is the 
total emission measure and the assumption / ~ 1 has been 
made. The volume of the emitting material is estimated from 
the measured flare area as V ~ A 3/2 . We apply Equation (10) 
to the eight well-constrained microflares, utilizing the best-fit 
Epstein DEM model in all cases. The total thermal energy U± is 
found to range between 3.9 x 10 28 — 1 .5 x 10 29 (see Table 2 for 
individual estimates). These estimates are similar to the mean 
value established by Hannah et al. (2008a) with an isothermal 
approach. The volume estimates obtained for these events are 
also comparable with Hannah et al. (2008a), ranging between 
2.4 x 10 27 -9.4 x 10 27 cm 3 . 

Comparison with the earlier derived values of L rad reveals 
an inconsistency. If a plasma is radiating simply, then the total 
radiated energy should not be less than the peak instantaneous 
thermal energy of the plasma. Indeed, during a recent study by 
Emslie et al. (2012) of 38 large solar eruptive events, it was 
noted that in general the total radiated energy from an event 
exceeded the peak thermal energy by a factor of ~3. 

For our measurements we find that L rad U±, with the 
thermal energy exceeding the total radiated energy by around 
an order of magnitude. This occurs despite our estimates of the 
emission measure and volume being broadly consistent with 
those found by Hannah et al. (2008a). 


The two main explanations for this energy inequality are 
that either thermal conduction losses are strong compared to 
radiative losses in the observed microflares, or that the filling 
factor / assumed to be unity for these estimates, is in fact small. 
A recent study of large flares by Guo et al. (2012) has suggested 
that 0.05 < / < 0.78, however, this may not hold for smaller 
events such as microflares (Moore et al. 1999). If we assume for 
our events the simplest relation between the thermal and radiated 
energy, assuming conductive losses to be small, that L ra d ~ U t h, 
then we can estimate lower limits on the filling factor /. For each 
event, we estimate / based on this assumption, finding a mean 
value of / ~ 5 x 10 -3 . The lower limits on /are similar to 
the estimates of the filling factor made by Baylor et al. (2011) 
in another recent microflare study, indicative of the idea that 
microflares and large flares may have different characteristic 
filling factors. In reality, however, the conductive losses in 
microflares are likely to be significant, particularly during the 
rise phase of the emission, accounting for at least some of 
the difference between L rad and U t h. To investigate this, we 
compared the contribution of non-thermal energy U n ±/ U t h (see 
Figure 9 and Table 2) in each event with I/ t h/T ra d- Weak non- 
thermal emission would suggest that more thermal conduction is 
taking place, as a mechanism is needed to evaporate hot material 
into the corona. Since this is not accounted for by L rad , one 
might expect a larger value of I/ t h/Trad in these cases. However, 
perhaps due to the small sample size, no relationship was found 
between these two measurements, though it must be recognized 
that the estimates of U n ± are only lower limits for each event. 

5.3. Non-thermal Energy 

A well-constrained model for the microflare thermal energy 
allows for a more accurate measurement of the non-thermal 
energy content by considering the remaining, high-energy, 
portion of the X-ray spectrum as non-thermal emission. Given 
a thick-target bremsstrahlung model, the power in non-thermal 
electrons above the low-energy cutoff is given by (e.g., Fletcher 
etal. 2013), 

8-1 F 

Pn±(> Ec) = y-^:— (ID 

0 — Z Ej C 

where P nt h is in erg s _1 , 8 is the power law spectral index of 
electrons, and F is the number of electrons per second above the 
low-energy cut-off E c . 

As discussed by many authors (see Holman et al. 2011), the 
low-energy cutoff E c can only be estimated, as it is masked by 
the thermal emission that dominates at lower energies. Hence, 
the values of E c used in Equation (11) should be consid- 
ered upper limits, which lead to lower limits on the amount 
of non-thermal power. The total non-thermal energy is esti- 
mated by integrating the power in non-thermal electrons over 
the non-thermal emission duration, and may also be considered 
lower limits. For a typical microflare, the non-thermal emission 
duration covers tens of seconds, while large flares may pro- 
duce non-thermal energy over several minutes. The values of 
non-thermal energy are obtained by integrating the estimated 
power in non-thermal electrons over the 60 s fitting duration. 

We find that the non-thermal power estimates range from a 
minimum of ~1.2 x 10 26 erg s -1 for event #10 to a maximum 
of ~1 .35 x 10 27 erg s -1 for event #2. The exception is event 
#5 which shows negligible non-thermal power when fitting the 
multi-thermal model. These estimated powers can be compared 
with the results of Hannah et al. (2008a), who estimated the 
power in non-thermal electrons for 4236 microflare events under 
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Figure 9. Thermal and non-thermal energy estimates for each microflare under (a) a multi-thermal assumption and (b) an isothermal assumption. 

(A color version of this figure is available in the online journal.) 


1.6 


an assumption of isothermal plasma, finding a mean power of 
P(> E c ) ~ 10 26 erg s -1 with a 2cr range of 10 25 -10 28 erg s _1 . 

Multiplying our estimated non-thermal powers by the 
RHESSI integration time, we find the estimated non-thermal 
energy content for these microflares is in the range 7 x 
10 27 -8 x 10 28 erg, again excluding event #5 where the estimated 
non-thermal power is negligible. The values for each event are 
listed in Table 2. 

5.4. Energy Balance and Comparison with 
an Isothermal Assumption 

For comparison, the thermal and non-thermal energy content 
for these microflares was also estimated following the approach 
of Hannah et al. (2008a), using an isothermal assumption with 
the addition of a thick target bremsstrahlung component to 
model non-thermal particle emission. These results are shown 
in the rightmost columns of Table 2. 

Here, we find that the estimated thermal energy content 
remains comparable with the estimates obtained under the multi- 
thermal model (see Figure 9). The similarity in thermal energy 
estimates may be explained by the distribution of contributions 
to the total instantaneous thermal energy in the multi-thermal 
case. Despite the multi-thermal nature of the plasma, most of the 
thermal energy is due to high-temperature plasma constrained by 
RHESSI observations. To illustrate this, in Figure 10 we integrate 
over the best-fit DEM obtained for microflare #1, gradually 
increasing the upper temperature limit of the integral. Here, the 
emission below log T ~ 6.5 contributes only 10% of the total 
thermal energy estimate for this event, while almost 50% of the 
total energy comes from emission above log T ~ 7.0. 

By contrast, the amount of non-thermal energy is substantially 
increased in an isothermal scenario compared with a multi- 
thermal model (Figure 9). The estimated power in non-thermal 
electrons in this model varies between 1.2 x 10 26 erg s -1 (event 
#9) and 3.2 x 10 27 erg s -1 (event #4). Although above the 
mean power estimated in Hannah et al. (2008a), these values 
are within their established 2 a range. 

For both the isothermal and the multi-thermal fits, we compute 
the fraction U n ±/U±. Figure 9 illustrates the change in energy 
balance that results under a multi-thermal assumption; the 
minimum non-thermal energy is in general approximately 
30% of the estimated thermal energy. By contrast, under 
a traditional isothermal assumption, the minimum estimated 
non-thermal energy is usually comparable with the thermal 
energy, the exception being event #9 where U nt h/U t h ~ 0.13. 
This change is due to the multi-thermal spectrum pushing 
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Figure 10. Contribution to the total instantaneous thermal energy as a function 
of maximum integration temperature for a multi-thermal DEM given by 
Equation (2). In this example, the best-fit DEM from the 2011 June 5 microflare 
(#1) was used. The curve in the figure is obtained by integrating over this best- 
fit DEM from low T, gradually increasing the upper limit of the integral, and 
comparing the obtained energies to the total energy obtained by integrating 
over all T. 

(A color version of this figure is available in the online journal.) 

the thermal component and thus the low-energy cutoff E c to 
higher energies (see Table 2). Here, the change in E c can be 
up to several keV, strongly affecting the estimated non-thermal 
energy contribution. It should be emphasized however that these 
are only lower limits on the level of non-thermal energy; the 
maximum cannot be constrained due to the uncertainty in E c . 

6. DISCUSSION AND CONCLUSIONS 

In this paper, we have, for the first time, combined EUV 
data from AIA and X-ray spectra from RHESSI to carry out 
joint forward- fitting of the DEM distribution for a selection of 
microflares. Although the sample size of events here is small, 
the results presented here are strong and have implications for 
future studies of flare DEM distributions and energetics. A study 
with a larger number of events is planned. We summarize our 
results below. 

First, we have found that a single Gaussian DEM distribution, 
although effective on AIA data in isolation, is unable to jointly 
fit AIA fluxes and RHESSI spectra for any of the studied 
microflares. This is due to the fact that such a function cannot 
decrease sufficiently steeply at high temperatures as required by 
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the RHESSI X-ray observations while simultaneously remaining 
broad enough to accurately reproduce AIA observed fluxes at 
low temperatures. We therefore conclude that a single Gaussian 
DEM distribution should not be used in microflares. 

Based on this finding, we tested the hypothesis that a uni- 
form DEM with a high-temperature cut-off is an adequate 
distribution to describe the emission measure distribution of 
microflares. This hypothesis was tested using an Epstein pro- 
file with steepness parameter n = 10 (see Equation (2)). It was 
found that, for 80% of the studied events, this simple model 
was able to fit both the AIA and RHESSI observations. Fur- 
thermore, for one event, #5, the combined EUV and X-ray 
emission could be fit using only a multi-thermal plasma, sug- 
gesting a very weak or non-existent non-thermal component 
(see Figure 5). 

We conclude that for most cases, 80%, a simple boxcar dis- 
tribution is a sufficient first-order approximation for describ- 
ing microflare DEM within the uncertainties of the AIA and 
RHESSI instruments. However, for two of the studied events, 
a joint fit to the EUV and X-ray data was not possible using 
either the Gaussian or Epstein DEM models. Consequently, the 
use of more complex DEM functions, such as a multiply peaked 
distribution, is justified in these cases. 

The thermal energy and the total radiated energy of the 
microflares was investigated. It was found that, when a filling 
factor of / ~ 1 was assumed, the thermal energy exceeded 
the total radiated energy by two orders of magnitude. This 
was corroborated by estimating the radiated energy using 
background-subtracted GOES X-ray data for each event (White 
et al. 2005; Emslie et al. 2012). We can postulate that either 
conductive losses are dominating the cooling of the coronal 
plasma for these events, or the filling factor is significantly less 
than unity. By neglecting conductive losses and postulating that 
Trad ~ Uth, we can estimate lower limits on the filling factor for 
these events to be / ~ 5 x 10 -3 . These estimates are plausible 
in view of the uncertainty surrounding the filling factor of the 
plasma (e.g., Cargill & Klimchuk 1997), although in reality 
some contribution to energy losses from conduction is expected 
during microflares, particularly during the rise phase of these 
events. Such low filling factors also appear to be inconsistent 
with more recent estimates of/ provided by Guo et al. (2012), 
who suggested that/ should be of order 0.1 or 1. However, Guo 
et al. (2012) studied only large events. Typical filling factors for 
microflares may be different from those for large flares (e.g., 
Moore et al. 1999; Baylor et al. 2011). 

The non-thermal energy was also investigated. In Section 5, it 
was shown that the application of a multi-thermal model with a 
thick-target non-thermal component substantially reduces the 
lower limit on the level of energy in non-thermal particles 
during microflares (Figure 9) compared to assumption an 
isothermal component. In general it was found that the minimum 
non-thermal energy content required was approximately 30% of 
the estimated thermal energy. This is in contrast to the traditional 
isothermal model, where the thermal and non-thermal energies 
were found to be comparable. With this model therefore there 
is much more thermal energy compared to non-thermal energy 
though it must be noted that the upper limit on the non-thermal 
energy remains unconstrained since the low-energy cutoff is not 
known. This suggests, though, that the standard model whereby 
the thermal emission in flares observed is due to chromospheric 
evaporation caused by the energy deposition of accelerated 
particles may not be appropriate in microflares. Instead, another 


source of energy, such as direct heating, accelerated protons, 
plasma waves, or DC-electric fields, may be the cause of the 
observed thermal energy. 

The authors are grateful to Dr Markus Aschwanden for 
helpful discussions and for making his software available within 
SSW/IDL, which provided the foundations of this work. We 
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flare energetics. 
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